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ABSTRACT 

A  numerical  method  of  calculating  deformations  and 
stresses  in  an  elastic  solid  propellant  grain  having  a  circular  perfora- 
tion  and  flat  ends  and  bonded  to  a  rigid  motor  case  was  used  to  calcu¬ 
late  stresses  and  deformations  caused  by  axial  acceleration  loading. 

Two  grain  end  situations  were  considered:  both  ends  free  and  one  end 
free  with  the  other  bonded  to  the  motor  case.  The  dimensionless  dis¬ 
placements  of  critical  parts  of  the  grain  were  calculated  as  functions  of 
the  inner-to-outer  radius  ratio  a,  the  length-to-diameter  ratio  and 
Poisson's  ratio  v.  The  dimensionless  radial  and  shear  stresses  at  the 
propellant-motor  case  interface  were  also  calculated  as  well  as  the 
percentage  of  total  load  carried  by  the  bonded  end  for  grains  bonded  on 
one  end. 

The  equations  of  motion  of  a  viscoelastic  solid  were 
reduced  to  a  single  Poisson  equation  by  assuming  that  the  displacements 
do  not  vary  with  the  coordinate  direction  along  which  the  acceleration 
is  applied.  A  general  solution  to  this  equation  was  obtained  for  hollow 
cylinders  of  infinite  length  having  generators  parallel  to  the  direction 
of  acceleration  and  transverse  cross  sections  with  p(p  ^  0)  axes  of 
symmetry.  The  solution  was  applied  to  solid  propellant  grains  of 
infinite  length  having  star-shaped  internal  perforations  to  determine 
the  stresses  and  deformations  caused  by  axial  acceleration  under  the 
conditions  of  zero  displacement  at  the  outer  boundary  and  zero  surface 
stresses  at  the  inner  boundary. 
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DEFORMATIONS  AND  STRESSES  IN  AXIALLY  ACCELERATED 
CASE-BONDED  SOLID  PROPELLANT  GRAINS  OF  FINITE  LENGTH 


Charles  H.  Parr 

< 


1 .  Introduction 

A  numerical  method  for  calculating  deformations  and 
stresses  in  an  elastic  solid  propellant  grain  of  finite  length  having  a 
circular  perforation  and  flat  ends  and  bonded  to  a  rigid  motor  case 
was  reported  previously*.  While  the  method  considered  internal 
pressure,  shrinkage,  and  axial  acceleration  loads,  the  calculations 
presented  were  concerned  primarily  with  shrinkage  loads.  Addi¬ 
tional  calculations  using  this  method  have  been  made  for  axial  acceler¬ 
ation  loads.  Two  end  conditions  were  considered:  two  free  ends,  and 
one  fixed  and  one  free  end.  As  pointed  out  previously,  it  is  possible 
to  present  results  in  dimensionless  form  in  terms  of  only  three 
parameters:  the  inner-to-outer  radius  ratio  a,  the  length-to-diameter 
ratio  \,  and  Poisson' s  ratio  v. 

2.  Correction  to  Previous  Report 

An  incorrect  statement  was  made  on  page  11  of 
Report  P-61-17  concerning  the  state  of  stress  at  the  midsection  of  a 
cylinder  of  finite  length  which  is  geometrically  symmetrical  about 
the  midsection  and  loaded  by  axial  acceleration.  It  was  stated  that  the 
stress  condition  in  the  midsection  plane  was  exactly  the  same  as  in  an 
infinitely  long  cylinder.  Actually,  at  this  plane  the  radial  displace¬ 
ment,  normal  stress,  and  normal  strain  are  zero  as  in  the  infinite 
cylinder  as  will  be  shown  below.  The  axial  displacement  and  shear 
stress  and  shear  strain  at  the  midsection  plane,  however,  cannot  be 
evaluated  from  considerations  of  the  midsection  alone  and  must  depend 
on  other  properties  of  the  cylinder,  specifically,  the  length. 


*Rohm  &  Haas  Company,  Quarterly  Progress  Report  on  Engineering 
Research,  No.  P-61-17,  June  1962. 
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Consider  the  accelerated  cylinder  in  Fig.  1.  If  the 
acceleration  is  considered  positive  as  shown,  a  negative  acceleration 
must  be  equivalent  to  an  acceleration  in  the  opposite  direction. 

Because  linear  analysis  is  used,  a  simple  reversal  of  the  signs  of 
the  stresses,  strains,  and  displacement  results  if  the  direction  of 
acceleration  is  changed  .  Consideration  of  this  point  leads  to  the  con¬ 
clusion  that  the  radial  displacement  along  the  geometric  line  of 
symmetry  m-m  must  be  zero.  If  the  radial  displacement  is  zero, 
the  radial  and  circumferential  strains  are  zero  along  this  line  by 
definition. 

To  demonstrate  that  the  axial  strain  is  zero  is  more 
difficult.  From  the  equation  for  radial  displacement,  Eq.  9  in  Report 
P-61-17  (neglecting  shrinkage  terms) ,  namely 

^  +  [  1  -  v]  «)>]  -  V  p  q  Z  ,  (9) 

there  results  from  setting  the  displacement  equal  to  zero* 

+  [  1  -  v]  <|)  «  .  (  34) 

The  first  and  second  derivativ';s  of  Eq.  34  are 


2  V  pT)  Z 
1  +  V 


(35) 


and 


9^«t>  .  ZvTl  Z 


3p^  ^  ap*  ■  ^  + 


(36) 


Addition  of  the  field  equations  given  by  Eqs.  7  and  8,  namely 

^  -  i  |i  +  ^  -0 

ap*  P  ^  aq*  ■ 


(7) 


'Equations  are  numbered  consecutively  with  Report  P-61-17. 
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and 

8p*  ~  P  ^  8ti*  ■  8t,* 

and  substitution  of  Eqs.  35  and  36  gives 


(8) 


^  =  V  ^  .  (37) 

8ti2 

Because  of  the  anti- symmetric  quality  of  the  radial  displacement  about 
m  -  m  for  the  axial  acceleration  loading,  it  can  also  be  concluded  that 
82u 

- ^  =0  (38) 

8ti* 


along  m  -  m,  i.e.,  there  is  an  inflection  point  in  the  radial  displace¬ 
ment  at  m  -  m.  Eqs.  37,  9,  and  38  when  combined  result  in 


8ti^  8ii* 


along  m  -  m.  Eqs.  7  and  8  now  yield  for  the  general  form  of  <j)  and 
along  the  line  m  -  m 


4*  s  A  +  B 

and 

=  C  p^  +  D  , 

which  can  be  evaluated  with  aid  of  the  boundary  conditions 

u=u  sO  atpsl 

p  q 

and 

a  s  a  rO  atpso 

pTi  p  I' 

to  give 

and 

vp  =  [  p*  -  0*1 


(40) 


(41) 


(42) 


(43) 
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where  t)q  is  the  axial  coordinate  at  m  -  nn.  By  substituting  Eq.  42  and 
43  into  the  expression  for  axial  strain,  namely 


€ 


[1±A 

p 


(44) 


the  axial  strain  can  be  found  to  be  zero  along  m  -  m.  Thus  the  three 
normal  strains,  and  hence  also  the  three  normal  stresses,  are  zero 
at  the  midsection  as  well  as  the  radial  displacement. 

Physical  considerations  similar  to  those  made  pre¬ 
viously  for  the  radial  displacement  will  establish  that  under  axial 
acceleration  not  only  the  radial  displacement  but  all  of  the  normal 
strains  and  stresses  as  well  are  anti- symmetric  functions  with  respect 
to  the  midsection,  while  the  axial  displacement  shear  stress,  and 
shear  strain  are  symmetric  functions.  Thus,  it  is  possible  to  deduce 
the  complete  solution  of  a  cylinder  which  is  geometrically  similar 
with  respect  to  the  midsection  by  considering  only  half  of  the 
cylinder  in  the  numerical  solution,  with  resulting  significant  savings 
in  time  and  cost. 


3.  Alternative  Method  for  Obtaining  Axial  Displacements 

In  Report  P-61-17,  the  axial  displacement  relation, 
Eq.  10,  was  obtained  by  integration  of  the  axial  strain  €  . 


u  s  — 


1  +  V 


dri 


+  r|6  +  +  f(  p) 


(  10) 


It  was  necessary  to  obtain  the  function  f(  p)  by  an  integration  of  shear 
strain  relations  when  the  axial  displacement  was  not  known  in  some 
radial  plane. 

An  alternative  method  of  obtaining  axial  deflections 
from  the  integration  of  the  shear  strain  relation  has  proven  more 
useful  than  Eq.  10.  Elimination  of  \  from  the  shear  strain  relation 


du 


3u 


= 


2[1  +  v1  d 


(45) 


and  the  radial  displacement  relation 

Up  :  p  ^4*  +  [l-'']4>^  -  vpT)Z  , 


(9) 
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4 


I 


results  in 


.  LkLuI  _  [i.v]  +  vpZ  (46) 

which,  upon  integrating  subject  to  the  condition  =  0  at  p  =  1, 
produces 


P  911  J 


dp 


(47) 


4.  Results  of  Calculations  for  Axial  Acceleratiori  Loading 
4.  1  Two  Ends  Free 

Calculations  were  made  to  obtai-.  stresses  ar.d  dis¬ 
placements  in  a  grain  having  both  ends  free  (  free-free  cylinder) 
subjected  to  axial  acceleration  loading  (  see  Fig.  1)  ,  Results  could 
be  presented  in  dimensionless  form  with  only  three  independent 

di 

parameters  -  the  radius  ratio  a  =  ^  ,  the  length-to-  diameter  ratio 

A.  s  ,  and  Poisson' s  ratio  v.  The  displacements  are  represented 
in  dimensionless  form  as 

u  E  u  E 

r  ,  z 

-  and  -  , 

b%g 

where  u  and  u  are  the  actual  radial  and  axial  displacemer.r.s 
respectively,  E  is  the  elastic  mod'olus  of  the  propellant,  y  is  the 
propellant  density,  and  g  is  the  axial  acceleration.  Or.  most  of  the 
figures  the  points  actually  calculated  are  marked  b/  appropriate 
symbols  so  that  the  interpolation  and  extrapolation  of  the  curves 
would  be  evident.  On  some  figures,  symbols  were  omitted  for 
clarity. 

The  dimensionless  deflections  A,  B,  and  C  defined 
in  Fig.  1  are  plotted  in  Figs.  2  through  7  as  functions  of  a  and  for 
Poisson' s  ratio  v  =0.5.  The  midsection  displacement  A  approaches 
the  displacement  of  an  infinite  cylinder  for  large  len.gth-to-diameter 
ratios  as  shown  in  Fig.  2.  Note  in  Fig.  7  that  the  deflection  C  is  not 
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results  in 


which,  upon  integrating  subject  to  the  condition  s  0  at  p  =  1, 
produces 


u 


T1 


9kjj  [  1  -  V  ] 
8ti  “  p 


dp 


(47) 


4.  Results  of  Calculations  for  Axial  Acceleration  Loading 
4.  1  Two  Ends  F ree 

Calculations  were  made  to  obtai’i  stresses  ar.d  dis¬ 
placements  in  a  grain  having  both  ends  free  (  free  -free  cvli’^der) 
subjected  to  axial  acceleration  loading  (  see  Fig.  1)  ,  Results  could 
be  presented  in  dimensionless  form  with  or.ly  three  ir.depeader.t 

parameters  -  the  radius  ratio  a  z  ^  ,  the  ler.gth-to-  dis.meter  ratio 

A.  :  '2^  ,  and  Poisson' s  ratio  v.  The  disp.lacemerts  are  represented 
in  dim.ensionless  form  as 

u  E  u  E 

r  j  z 

-  and  -  , 

bSg  bSg 

where  u  and  u  are  the  actual  radial  ar.d  axi-il.  displacemer '.s 
respectively,  E  is  the  elastic  modulus  of  the  propellant,  is  the 
propellant  density,  and  g  is  the  axial  acceleratic”'..  Or.  most  of  the 
figures  the  points  actually  calculated  are  marked  V>/  appropr'aie 
symbols  so  that  the  interpolation  and  extrapolalio of  the  c'.rves 
would  be  evident.  On  some  figures,  symbols  were  omitted  for 
clarity. 

The  dimensionless  deflections  A,  B,  ar’d  C  defired 
in  Fig,  1  are  plotted  in  Figs.  Z  through  7  as  functions  of  a  and  \  for 
Poisson' s  ratio  v  =0.5.  The  midsection  displacement  A  approaches 
the  displacement  of  an  infinite  cylinder  for  large  length-to-diameter 
ratios  as  shown  in  Fig.  2.  Note  in  Fig.  7  that  the  deflection  C  is  not 
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a  monotonic  function  of  a.  The  variation  of  A,  B,  and  C  with  Poisson's 
ratio  for  particular  values  of  a  and  \  is  shown  in  Figs.  8,  9>  and  10. 

Figures  11  through  14  and  Figs.  15  through  18  illustrate 
the  effect  of  \  and  a  on  the  radial  and  shear  stress  respectively  at  the 
propellant  motor  case  interface  for  a  Poisson' s  ratio  of  v  :  0.  5.  A 
singularity  in  the  stress  distribution  occurs  at  the  corner  formed  by 
a  free  end  and  the  motor  case.  Thus,  any  stresses  calculated  near 
these  points  are  physically  meaningless  because  any  real  material 
will  deviate  from  the  assumed  model.  In  the  figures  illustrating  stress 
distribution,  the  stresses  near  this  singularity  are  indicated  by  dashed 
lines.  The  effect  of  this  singularity  on  the  numerical  solution  has  not 
been  fully  evaluated  but  is  believed  to  be  small  even  at  small  dis¬ 
tances  from  the  singular  point.  The  effects  of  Poisson' s  ratio  on  the 
displacements  near  the  end  are  complex,  but  at  points  removed  from 
the  ends  the  stresses  are  not  highly  dependent  on  Poisson' s  ratio. 

4.  Z  One  End  Fixed 

For  cylinders  having  one  end  bonded  ( fixed- free 
cylinder) ,  as  illustrated  in  Fig.  19,  results  similar  to  those  obtained 
for  the  free-free  cylinders  are  presented  in  Figs,  20  through  43.  Note 
that  the  proper  interpretation  of  the  sign  of  the  deflection  and  stresses 
depends  on  the  direction  of  acceleration  while  the  magnitude  is 
independent  of  the  direction  of  acceleration.  Figs.  20  through  25 
demonstrate  the  variation  of  the  quantities  A,  B,  and  C  as  defined  in 
Fig.  19  for  various  valuei^  of  \  and  a  but  only  for  Poisson' s  ratio  of 
0.5.  Figs.  26  through  27  demonstrate  the  variation  of  A,  B,  and  C 
with  Poisson' s  ratio.  Although  these  figures  are  not  a  complete 
parameter  study,  they  do  indicate  the  complex  nature  of  the  depend¬ 
ence  on  Poisson' s  ratio. 

Figs.  29  through  33  show  the  load  carried  by  the  end 
bond  as  the  percentage  of  the  total  load  carried  by  the  motor  case  in 
accelerating  the  grain.  Figs.  29  and  30  are  crossplots  in  a  and  A.  of 
this  data  for  v  s  0.  5.  Figs.  31  through  33  illustrate  the  dependence 
on  Poisson' s  ratio. 
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Fig.  I  Cross  section  of  cylinder  with  free  ends  deformed  by  acceleration. 
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Flg.  2  Axial  displacemeni  at  midseciion  of  free -free  cylinder,  v«0.5. 
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Axial  displacement  of  corner  of  free-free  cylinder,  i/=05. 
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Fig.  7  Radial  displacement  of  corner  of  free-free  cylinder,  (/‘OlS 
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Fig.  0  Maximum  axial  displacement  at  midsection  of  free -free  cylinder. 
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Fig.  9  Axial  displacement  of  corner  of  free-free  cylinder. 
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Rg.  10  Radiol  displacement  of  corner  of  free -free  cylinder. 
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z/b,  DISTANCE  FROM 
'  GRAIN  END 


Fig.  II  Dimensionless  rodial  stress  at  propellant- motor  case  interface,  \-0.5, 
v-0.5. 
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Rg.  12  Dimensionless  radial  stress  q+  propellant -motor  case  interface,  X=  1.0, 
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Eig.  14  Dimensionless  radial  stress  af  propellant -motor  case  interface,  X=4.0, 
V  *0.5. 
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Eig.  14  Dimensionless  radial  stress  at  propellant -motor  cose  interface,  X*4.0, 
i/»0. 5. 
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0  0.2  Q4  0.5 

z/b,  DISTANCE 
FROM  GRAIN  END 

Fig.  15  Dimensionless  shear  stress  at  propellont-motor  cose  interface,  X=0.5, 
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Fig.  16  Dimensionless  shear  stress  at  propel  lent -motor  case  interface  X  =  IO 

v'0. 5. 
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rig,  19  Cross  sections  of  cylinders  with  one  fixed  end ,  deformed  by  occeierotion. 
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Fig.  24  Radial  displacement  of  free  corner  of  fixed-free  cylinder,  i'  =  0.5. 
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25  Radial  displacement  of  free  corner  of  fixed-free  cylinder,  v=0.5. 
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r/b,  RADIUS 


Fig.  35  Distribution  of  shear  stress  at  head  end  bond,  w=0.5. 
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z/b,  DISTANCE  FROM  FIXED  END 

Fig.  36  Dimensionless  radial  stress  at  propellant -motor  case  interface,  X=Q5, 
i/=0.5. 
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Rg.  39  Dimensionless  radial  stress  at  propellant -motor  case  Interface,  X=4, 

V  =  0. 5  . 
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z/b  DISTANCE  FROM 
FIXED  END 


Fig.  40  Dimcnsiontess  shear  stress  at  propellant-motor  case  interface,  X=Q5, 

v-0.5. 
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42  Dimensionless  shear  stress  at  propellant -motor  case  interface,  \-Z.O.  v-O.b. 
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DEFORMATION  OF  VISCOELASTIC  CYLINDERS  OF 
INFINITE  LENGTH  AND  ARBITRARY  CROSS  SECTION 
UNDER  AXIAL  ACCELERATION  LOADS 


Charles  H.  Parr 


1.  Introduction 

A  knowledge  of  the  stresses  and  deformations  in 
solid  propellant  rocket  motors  due  to  axial  acceleration  loads  is 
necessary  for  analysis  of  motor  structural  integrity.  An  analysis 
of  the  effects  of  finite  length  on  the  stresses  and  deformations  in  pro¬ 
pellant  grains  having  circular  perforations  under  axial  acceleration 
is  presented  in  another  section  of  this  report.  This  section  deals 
with  the  axial  acceleration  of  propellant  grains  of  infinite  length  whose 
internal  perforations  are  not  circular  but  have  a  number  of  axes  of 
symmetry,  such  as  the  common  types  of  star  perforations. 

The  equations  of  motion  of  a  viscoelastic  solid  were 
reduced  to  an  uncoupled  set  by  assuming  that  the  deformations  do  not 
vary  with  the  coordinate  direction  along  which  the  acceleration  is 
applied.  These  results  were  then  applied  by  the  use  of  conformal 
mapping  to  axially  accelerated  cylinders  with  cross  sections  having  p 
axes  of  symmetry  ( p  ft  0)  .  A  closed  form  solution  was  obtained  for 
the  axial  displacement  and  then  evaluated  for  a  particular  set  of 
boundary  conditions  and  a  particular  geometry. 


2.  Equations  of  Motion 


The  equations  of  motion  for  a  viscoelastic  solid  may 

be  shown  to  be* 


+  G  [t"  u(t)]  +  X  s  Y  ^ 


(1) 


‘This  formulation,  with  the  exception  of  acceleration  and  body  force 
terms,  has  been  given  by  Elder  [  1]  . 
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+  Y  =  V 


d^v 

a 


and 


(2) 


K  + 


9  e(  t) 
dz 


+  G  [y*w(t)]  +  z  =  V  —  . 

8 


where  the  functional  notation  of  Volterra  [  2] 

G[f(t)]s  f(t)  G(0)  ^  [G(t-t,)]dti 


(3) 


(4) 


is  used.  Here 

G(  t)  ■  viscoelastic  shear  relaxation  modulus, 

K{  t)  ■  viscoelastic  bulk  relaxation  modulus, 

u,  V,  w  ■  displacements  in  direction  of  Cartesian  coordinates, 
X,  y,  z  respectively*  , 

y  z  material  density, 

t  ■  time. 


8u  ,  8v  ,  8w  ,,,  , 

e  :  fT-  +  s—  +  T"  ~  dilatation,  and 
8x  8y  8z 

X,  Y,  Z  :  body  forces  in  the  x,  y,  z,  directions,  respectively. 

It  should  be  noted  that 

G[v*u(t)]  :  V"G[u{t)]  .  (5) 

Under  the  restriction  that  u,  v,  and  w  are  invariant 
with  z,  Eqs,  1,  2,  and  3  reduce  to 


*The  displacements  u,  v,  and  w  are,  in  general,  functions  of  x,  y,  and 
z  as  well  as  t.  For  brevity,  the  dependence  on  the  coordinates  is  not 
expressed  in  function  form  but  the  dependence  on  t  is  sometimes 
expressed  to  aid  in  expressing  the  Volterra  and  Laplace  functionals. 


-52- 


K  + 


9^u(  t)  d^v(  t) 
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+  X  =  V 


a^u 
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(6) 


K  + 


t)  d^v(  t) 

9xay  2 

ay 


+  G|  t) 


ax' 


ay' 


a^v 


+  Y  =  y  ^  .  (7) 

at^ 


and 


a^w(  t)  ^  a^w(  t) 


ax' 


ay' 


+  Z  -  Y 


a^w 

at^ 


(8) 


Eqs.  6  ar.d  7  are  now  lincoupled  from  Eq.  8.  Eqs. 

6  and  7  correspond  to  the  usual  plane  strain  equations  of  visco¬ 
elasticity  and  may  be  solved  separately  for  any  given  boundary  condi¬ 
tions  and  body  forces.  Henceforth,  consideration  will  be  given  to  the 
solution  of  Eq.  8  with  geometries  limited  to  cylinders  having  genera¬ 
tors  parallel  to  the  z-axis.  Solutions  fcr  the  coupled  equations,  Eqs, 

6  and  7,  will  not  be  considered  in  this  analysis  but  can  be  superposed 
with  solutions  of  Eq.  8.  Note,  however,  that  for  a  loading  of  only 
acceleration  in  the  z-direction,  the  solution  of  Eqs.  6  and  7  is  trivial, 
i.  e, ,  u  a  V  :  0. 

Body  forces  will  be  restricted  to  the  weight  of  the 
body.  It  is  evident  that  the  weight  per  unit  volume  can  be  expressed 
as  the  product  of  the  density  y  a  pseudo  acceleration,  the 
gravitional  constant,  and  can  thus  be  included  in  the  right  side  of 
Eqs.  6,  7,  and  8.  Consequently  the  body  force  will  no  longer  be 
explicitly  considered. 

Boundary  conditions  applicable  to  Eq.  8  may  consist 
of  the  specification  of  the  displacement  or  shear  stress  as  prescribed 
functions  of  time.  Typically  these  may  take  the  form 

w(  t)  =  f{  t)  on  B,  (  9) 


Tgz(t)  =  h(t)  on  B^ 


(  10) 


and 


-53- 


where  Bj  and  B2  are  the  boundaries  of  a  hollow  cylindrical  body  and 
s  denotes  the  normal  to  the  surface. 

The  displacement  w  can  be  considered  to  consist  of 
two  parts,  Wj  and  Wj ,  where  wj,  a  function  of  time  only,  is  the  dis¬ 
placement  of  boundary  Bj  and  w^  is  the  displacement  relative  to  the 
boundary  Bj .  The  displacement  Wj  may  be  associated  with  a  rigid 
body  displacement.  The  displacement  W2  is  associated  with  the  defor¬ 
mation  of  the  cylinder  material  and  is  a  function  of  the  space 
coordinates  x  and  y  and  of  time  t.  The  displacement  W2  will  be  further 
restricted  by  neglecting  dynamic  effects  so  that  the  magnitude  of  the 
variation  of  W2  with  time  is  much  less  than  that  of  Wj  .  Thus 

9^w  ^  9^w , 

9t^  9t^ 

and  since 


9w  1  9w 
8x  ■  9y 


Eq.  (  8)  may  be  written 


+  t) 

9y" 


V 


9^Wi(  t) 

9t2 


(11) 


Specification  of  the  overall  body  acceleration  may 
now  be  made  independently  of  the  displacements  associated  with 
deformations. 


To  obtain  a  solution  to  Eq.  11,  it  is  convenient  to  use 
the  Laplace  transform  with  respect  to  time  to  reduce  the  problem  to 
an  associated  elastic  problem.  Applying  the  Laplace  transform  to  the 
functional  G  r  W2(  t)  1  there  results 


:  G(  0)  W2{  p)  +  W2(  p)  p  G{  p)  -  G(  0)] 


=  P  W2(  P)  C(  p)  ,  ( 12) 

where  p  is  the  transform  variable,  not  to  be  confused  with  the 
number  p  of  axes  of  symmetry  of  the  transverse  cross  section. 
Expressing  the  acceleration  of  the  body  as  a  specified  function  A(  t) , 
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:  A(t)  .  (13) 

at* 

where 

A{  t)  =  0,  t  <  0  , 

Eq,  11  becomes 

y*  p)  :  ( 14) 

P'J^(p) 

or 

V*W2(p)  :  YP^(p)  J(p)  (15) 

since 

pC(p)  =  — 

P  J(  P) 


where  J(  p)  is  the  Laplace  transform  of  the  shear  compliance. 

Taking  the  inverse  Laplace  transform,  there  results 

y*L'‘{wj(p)^  sL‘‘^YP^(p)  J(p)^  *  (16) 


where  L 


{ 


^  denotes  the  inverse  Laplace  transform. 


Letting 


L'*(w^(  p)} 

L'*  {  YP7(P)  ^(P)} 


(17) 


Eq.  16  becomes 


^  ^  I 

ax*  ay* 


(18) 


Note  that  ijj  is  a  function  of  x  and  y  only  and  is  not  a  function  of  time. 
Eq.  18  is  valid  regardless  of  the  time  dependence  of  the  acceleration 
and  regardless  of  the  time  dependence  of  the  shear  compliance. 
Further  it  can  be  seen  from  Eq.  18  that  the  displacement  function  il* 
defined  by  Eq.  17  depends  solely  on  the  space  coordinates  x,  y  and  the 
boundary  conditions. 
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If  the  acceleration  A(  t)  is  applied  as  a  step  function, 
A(t)  =  gH(t)  , 

where  H(  t)  is  the  unit  step  function  and  g  is  a  constant,  then 
^(p)  •  g/P  .  and 


L'^(w2(  P)}  W2 
■‘(gvKp)}  ■ 


(19) 


Thus  for  a  step  function  in  the  acceleration,  the  axial  displacement 
at  each  point  is  directly  proportional  to  the  creep  compliance. 

3.  Application  to  Star  Geometries 

The  solution  of  Poisson' s  equation,  Eq.  18,  was  con¬ 
sidered  for  star  geometries  which  are  applicable  to  solid  propellant 
grains.  Consider  the  shape  shown  in  Fig.  1.  The  internal  perforation 
consists  of  p  branches  or  star  points  and  the  external  boundary  is 
circular.  Wilson  [  3,  4]  demonstrated  the  mapping  of  such  regions  by 


Rg.  I  Typical  star  perforated  propellant  grain. 
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fitting  the  internal  perforation  with  a  mapping  transform  which  maps 
the  star  shape  in  the  z*  or  real  plane*  onto  the  unit  circle  in  the  ^  or 
transformed  plane.  Likewise,  the  circle  defined  by  p  s  F  in  the  ^ 
plane  corresponds  to  an  irregular  line  in  the  z*  plane.  However,  if  F 
is  sufficiently  large  the  corresponding  contour  in  the  z’*‘-plane  is 
sufficiently  regular  to  be  considered  circular. 

The  mapping  function  used  by  Wilson  is  given  by 


=  u{  4)  =  A  4  + 


B 

rP- 


D 

.5p-l 


(20) 


where  p  is  the  number  of  axes  of  symmetry  of  the  grain  cross  section. 
Kantorovich  and  Krylov  [  5]  show  that  under  such  a  transformation, 
the  Laplace  operator  transforms  into 


/  3*  ,  3*  1  .  1  I 

a* 

+  -^1 

lax*  ax*  1  1“ '  ( 4)|  ^  1 

ag* 

ari"  / 

where  «'  (  4)  is  the  derivative  of  to  with 
21  in  polar  coordinates  by  means  of 

respect  to  4.  Expres 

*  .  .  i0 

4:§+iT)spe  , 

the  Laplace  operator  becomes 

/  a"  ,  aM  _  1 

1  ^ 

+  1  ±  +  ±  ill 

1  ax*  ay*  1  |w'  ( * 

1  ap* 

p  ^P  p*  ae*/ 

(21) 


Using  Eq.  2Z,  Eq.  18  may  be  written 

^  +  i  +1  ^ 

ap*  P  ^  p*  30*  ' 


(23) 


Now  consider  the  general  form  of  w(  4), 


*To  avoid  confusion  with  the  z  coordinate,  an  asterisk  is  used  to 
denote  the  complex  variable  z*  :  x  +  iy  . 
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<4i)  =  • 

Then 

«'(&)  =  ^(i-np) 


(24) 


or 


w'(  4) 


( l-np)  C 


~  ^  cos  np0  -  i  sin  npsj  . 


Then 

l«'(4)r 


(l-np)  C 

-  cos  npe 

np 

P 


(25) 


(l-np) 


-l2 


sin  np6 


(26) 


or,  upon  expanding, 


I «'(;)!  *  = 


(l-np)(l-mp) 


(  m+n)  p -  mp0  cos  np0 


•tt 


( l-np)  ( 1-mp)  C  C 

_ _ _  sin  mp0  sin  np0  , 

(  m+n)p 


(27) 


which  simplifies  to 
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2  _ 


( 1-np)  (1-mp)  C  C 

- (n^n)r-  -  P« 


(28) 


Let  ksm-n.  Then,  after  considerable  manipulation,  Eq.  28  may  be 
rewritten  in  the  form 


(l-np)*  C* 

.2np 


+  2 


cos  kp8 


( 1-np)  ( 1-np-kp) 
„(  2n+k)p 


(29) 


The  solutions  to  Laplace' s  equation  which  will  form 
the  complementary  part  of  the  solution  to  Eq.  23  are  well  known^ 
Those  forms  suitable  for  the  problem  under  discussion  are 


( +  Lj^  p  )  cos  kp0  +  U  in  p  +  V  , 


(30) 


where  ,  Lj^,  U,  and  V  are  to  be  determined  from  the  boundary 
conditions. 

The  particular  solution  of  Eq.  23  may  be  found  by 


considering 

s  p“  cos  kpO  . 


(31) 


Substituting  this  relation  together  with  Eq.  29  into  Eq.  23  and  equat¬ 
ing  coefficients  of  cos  kp6  results  in 


*See,  for  instance.  Moon  and  Spencer  [  6]  ,  p.  14. 
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a  s  2  -  (  2n  +  k)  p  , 


M 


C  C 
n  n+k 


nk 


and 


k  ^  0 


(32) 


The  particular  solution  is  then  formed  by  a  double  series  in  n  and  k 
of  terms  given  by  Eq.  31. 

The  complete  solution  is  then  given  by  adding  the  com¬ 
plementary  and  particular  solutions  which  results  in 


p2(l-np)  +  uin  p  +  V 


^n  ^n+k  2-(2n+k)pl  . 

- 5 -  p  '  '  ^  )  cos  kpQ  .  (33) 


This  is  the  solution  for  the  axial  displacement  function  in  an  infinite 
cylinder  whose  contour  in  the  transverse  cross  section  can  be  defined 
by  a  mapping  function  of  the  form  Eq.  20  (  or  Eq.  24)  . 

4.  Satisfaction  of  Boundary  Conditions 

Consider  the  boundary  conditions 

w(x,  y,  t)  =  wi(t)  +  W2(x,  y,  t)  =  f(  t)  on  B,  (9) 

and 

’’sz  '  on  ,  ( 10) 


By  definition 

w,(Bi,t)  =  0 


(34) 
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so  that  the  first  boundary  condition,  Eq.  9,  can  be  written 

w(x,  y,  t)  =  W|{t)  s  f(t)  =  A(  t)  dtdt  on  .  (35) 

Considering  the  second  boundary  condition,  Eq.  10,  where  s  is  the 
outward  normal,  the  shear  stress  can  be  expressed  in  terms  of 
shear  strain  and,  in  turn,  in  terms  of  the  axial  displacement  Wj  and 
normal  displacement  u  resulting  in 


T  r  G 
sz 


9u  ( t)  a  / 

,  OwJt) 


-  n 

dW2(  t) 

L  J 

9z  '  5s 

Taking  h(  t)  :  0,  boundary  condition  Eq.  10  thus  reduces  to 


(36) 


9w2(  t) 
ds 


=  0 


on  B, 


(37) 


Eq.  37  is  a  nomogeneous  Volterra  integral  equation  and  the  only 
continuous  solution  is  zero  [  7]  .  Thus 

=  «  • 

For  the  mapping  function  used,  it  can  be  shown  that 


_9_ 

ds 


1 


|w'  (  U  I  9p 

on  an  internal  boundary  and 


ds 


1  _d_ 

|w''  {  V  1  9p 


(38) 


(39) 


(40) 


on  an  external  boundary. 

Defining  B|  as  the  outer  boundary  and  B2  as  the 
inner  boundary,  the  boundary  conditions  in  the  transformed  plane, 
with  the  use  of  Eqs.  34,  38,  and  39,  become 


and 


W2(  t)  :  0  on  B  j 


0  o„B, 


(41) 


(42) 
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Since  these  conditions  are  independent  of  time  they  may  be  applied 
directly  to  the  function  i)j  to  yield 

41  s  0  on  Bj  ( 43) 

and 

^  ;  0  on  B2  .  ( 44) 

Letting 

p  ;  1  on  Bj 

and 

p  r  /3  on  B2 

the  unknowns  ,  L^^ ,  U,  and  V  in  Eq.  33  can  be  evaluated  with  the 
aid  of  Eqs.  43  and  44  allowing  Eq.  33  to  be  written  as 


j  2-np-kp  ^  ^2(  l-np)j  I^Zkp  _  pZkpj 


i  +  ;3 


ikp 


p~^^  cos  kp0  .  (  45) 


5.  Evaluation  of  Displacements  and  Shear  Stresses 

Due  to  the  form  of  the  mapping  function  co(  0  ,  Eq.  24, 

i9 

an  inverse  is  not  easily  obtainable.  Thus  the  value  of  4  =  p  e 
corresponding  to  a  given  location  z*  =  x  +  iy  cannot  be  explicitly 
obtained  but  must  be  approximated.  A  straightforward  way  of  evalu¬ 
ating  the  expression  for  displacement  is  to  perform  the  evaluation 
of  vl'  and  z*  for  a  series  of  p' s  and  0' s  and  then  use  interpolation  to 
obtain  values  of  displacement  at  specific  values  of  z.  In  this  manner 
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the  dimensionless  displacement  t{'  be  evaluated  at  any  point  for  a 
given  geometry  from  Eq»  45  by  specifying  the  and  /3.  The  actual 
time -dependent  displacement  can  then  be  obtained  by  performing 
the  operations  indicated  by  Eq,  17  once  the  acceleration  A(  t)  and  the 
creep  compliance  J(  t)  are  specified. 

The  shear  stress  is  somewhat  more  difficult  to 
evaluate.  Since  the  displacements  u  and  v  in  the  x  -  y  plane  are  zero 
for  a  loading  of  axial  acceleration,  the  shear  stresses  due  to  the 
axial  acceleration  load  are  given  by 


T 


iz 


(46) 


where  i  is  any  desired  line  lying  in  the  or  x  -  y  plane.  Taking  the 
Laplace  transform  of  this  expression  with  respect  to  time, 

=  pG(p)  ,  (47) 

the  transformed  shear  stress  can  be  written,  with  the  use  of  Eq.  17, 
as 

T^Jp)  :  Y  S(P)  If 
or  simply 

=  VA(.)|^  .  (48) 


Usually  the  maximum  v^ai  j.es  of  shear  stress  at  a 
point  are  wanted.  These  may  be  obtained  in  the  following  manner. 
The  relation  of  the  derivative  along  a  line  i  in  z^-plane  to  the 
derivative  along  a  corresponding  line  A  in  the  ^-plane  is  given  by 

^  Ia  • 

In  terms  of  the  p-0  coordinates 


(50) 
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But 

=  cos  a  (  51) 

and 

=  sin  a  ,  (52) 

where  a  is  the  angle  between  \  and  the  p  coordinate  line.  Substituting 
Eq.  51  and  £q.  52  in  Eq.  50  and  maximizing  with  respect  to  a,  there 
results 


Eqs.  48,  49,  and  53  yield 
=  V  A(  tp 


max 


The  displacement  function  i)*  has  been  evaluated  for 
the  six-point,  star-perforated  shape  shown  in  Fig.  2  which  has  been 
mapped  with  a  50  term  mapping  function.  The  original  ballistic  con¬ 
figuration  is  indicated  by  the  dashed  lines.  The  curvilinear  coordi¬ 
nates  are  shown  over  a  part  of  the  cross  section,  demonstrating  the 
slight  mismatch  between  the  true  circular  outer  boundary  and  the 
curvilinear  coordinate.  A  contour  map  of  the  displacement  is  shown 
in  Fig.  3  for  a  step-function  acceleration  load,  i.  e. ,  A{  t)  s  g  H{  t)  . 
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